The Power of Various Real- Valued Quantum 

Queries 



The computation of combinatorial and numerical problems on quan- 
tum computers is often much faster than on a classical computer in num- 
bers of queries. A query is a procedure by which the quantum computer 
gains information about the specific problem. 

Different query definitions were given and our aim is to review them 
and to show that these definitions are not equivalent. To achieve this 
result we will study the simulation and approximation of one query type 
by another. While approximation is "easy" in one direction, we will show 
that it is "hard" in the other direction by a lower bound for the numbers 
of queries needed in the simulation. The main tool in this lower bound 
proof is a relationship between quantum algorithms and trigonometric 
polynomials that we will establish. 

1 Introduction to Quantum Queries 

Since Grover's discovery of the quantum search algorithm 0, the notion of 
query complexity has played an important role in quantum computation (for 
a thorough introduction to quantum computation see Nielsen, Chuang, 
Grover studied the problem to find an index j e {0, 1, . . . , — 1} so that for a 
function 



/(j) = 1 holds. He constructed a quantum algorithm that needs 0{'\/N) eval- 
uations of /, if exactly one index j with /(j) = 1 exists, vifhile an algorithm on 
a classical computer would need Q{N) evaluations. His quantum evaluation of 
/ was given by the transformation 
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This query can easily be seen to be equivalent to the query 

\j) \b) ^ \j) \b ® f{j)) = \j) \b + f{j) mod 2) (1) 

which was used by Beals et al. [2] and Nayak and Wu (7). For this query they 
proved important lower bounds for many problems that involve Boolean queries. 

All these queries refer to the following model of computation. Let Q/ denote 
a query transformation for a fixed function /, let Uk denote arbitrary unitary 
transformations not depending on /, and let lip) be a starting state. Consider 
the quantum algorithm 

UnQfUn-l...UiQfUo\i^). (2) 

It returns a resulting state \xf) after measurement with probability Pxf We 
identify this state \xf) with its index Xf and map it to the solution space by a 
computable mapping ip not depending on / on a classical computer. 

We briefly recall some definitions from information-based complexity (see 
the introductions and surveys Traub, Werschulz 13 , and Traub, Wasilkowski, 
Wozniakowski ^12,). These were transferred to quantum computers by Heinrich, 

whose approach we follow. Let the correct solution to every problem /, / 
from the problem class F, be given by a mapping S : F ^ G. If the algorithm 
returns a ^fiixf) G G for every input f £ F with probability p^j , so that 

3 

2^ Pxf > 4' 

X/:||S(/)-./>(x/)||<£ 

it is said to solve the problem for the problem class F. The best precision e that 
can be achieved by a quantum algorithm, which uses only n queries, is called 
the n-th minimal quantum query error e'^{S, F) . Conversely the quantum query 
complexity comp^(5', F) is the minimal number of queries that have to be used 
in a quantum algorithm to get the desired result with precision e. 

With this formal framework we are able to deal with numerical problems, 
if we define what the query Q f is in this case. Since numerical queries have 
more than just two possible outcomes we have to extend our notion of a query. 
We will follow the formal approach of Heinrich to discuss the different queries. 
Suppose we want to query a function f : D K. To ease our discussion let us 
restrict to K C [0, 1]. We have to choose a mapping r : {0, 1, . . . , 2" — 1} — > D, 
decoding the input, and a mapping (3 encoding the output. 

There are two different approaches to encode a value in [0, 1] in a quantum 
computer. The first approach was established by Abrams and Williams |T], and 
Novak [I^, who proposed the following query 

\j) |o> = v/i-W(t(j)))Ij) |0) + V/3(/(t(j))) \j) |1) 

, , (3) 

gph-|j-) |i) = -^mAMlj) |0> + v/i-/3(/(r(j)))b-> |1> 

This query continuously changes between the two states |0) and We will 
refer to this query as the phase query since it can be viewed as a (relative) phase 
rotation in the space spanned by |0), In ^ the mapping /? was chosen as 

/3:[0,1] -> [0,1] 
X 1-^ x'^ 
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while in [ini /9 = id. 

Another approach to extend the Boolean query was chosen by Heinrich, [Sj . 
He encoded the values in [0, 1] via /3 to {0, 1, . . . , 2"^ — 1} and defined: 

Qf \j) \x) = \j) \x ® /3(/(r(j)))) = \j) \x + PifiriM mod 2™) (4) 

For the rest of this paper we will refer to this query definition as the bit query. 

Consider the problem to compute the mean. Let F := {/ : {0, . . . , iV — 1} ^ 
[0, 1]}. The solution operator Sn : F ^ [0, 1] for this problem is 



Brassard et al. pi, Novak ^OIj and Heinrich showed that for this problem 

comp?(5Ar,F) = e(i) 

with phase queries as well as with bit queries, where in both cases the upper 
bounds were shown by the amplitude estimation algorithm. 
Does this equivalence hold for all problems? 



2 Bit queries can simulate phase queries 

Suppose we are able to construct a unitary transformation that realizes a phase 
query by the use of bit queries. If we additionally show that the required number 
of bit queries is a small constant, no algorithm that uses phase queries could be 
(asymptotically) faster than an algorithm with bit queries. 

Definition 1. If it is possible to find unitary transformations Ui not depending 
on f with either 

Qf ^Un,QfUn,-i...UxQfUo (simulation) (5) 



or 



Qf — Un^QfUjig-i ■ ■ -UiQ^Uq <d (approximation) (6) 

for all f G F , we say that the query Qf is simulated or approximated up to S 
by Qf using Uq queries Qf . 

What does this definition of approximation say about the probability to 
measure the correct state? As it is shown e.g. in |H], box 4.1, suppose the 
probability that on input the state \lp) is measured with probability pu^ if 
U is applied, and with pv , if V is applied. Then 

\pu-Pv\<m-V\\. (7) 

Proposition 2. Two applications of Q^^"^ suffice to approximate Q^f"'^'^ by Q^*. 
The precision of the approximation only depends on the encodings /J^'"***^, Z?****, 
and the number of bits in the bit query Q^j^*" . 
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Fix the encodings f3^'^ : [0, 1] ^ {0, . . . , 2™ - 1}, and (3^^^^° : [0, 1] [0, 1]. 
Define : {0, . . . , 2™ - 1} -> [0, 1] as a mapping that satisfies /3'''* o /3li't = id. 
Let the starting state 1-0) be 

1^) = E io> + E ii) e ® 

We just state the algorithm for the basis states \j) |0), and |j) |1), the extension 
to arbitrary states follows from linearity. 

• Starting state is \j) |0) or \j) |1) . 

• Append |0) |0) G i/„ «) H.^: 

|j)|0)|0)|0) or |j)|l)|0)|0) 

• Copy (by addition modulo 2") the first register to the third register: 

|j)|0)b-)|0) or |j)|l)b-)|0) 

• Apply to the third and fourth register: 

|j>|0)b-)|/3^"(/Wj)))) or b-)|l)|j)|/3'^"(/(r(j)))) (8) 

• Define a mapping U by 

U \j) |0) \j) \x) =^l-/3phasc(^bit(^)) 1^-^ 

+ J/3Phase(^bit(^)) 1^-^ 1^-^ 

\ (9) 

U \j) |1) \j) \X) = - y'/3Phase(^bit(^)) 1^-^ 1^-^ 
+ y^l-/?Pha-(/3bit(^))|^-) |1) 

for every j e {0, . . . , 2" - 1} and x G {0, . . . , 2™ - 1}. Note that U does 
not depend on /. Apply U to the state resulting in: 

^l_^phasc(^bit(^bit(J(^(_^-)))))|j-) |0) 1^-) |^bit(^(^(^.)))^ 

+ ^f3P^^^{l3^j\f3^-^Hf{Tm)) \J) |1> |J> (10) 

or 

- y^/3ph-c(^bit(^bit(j(^(^-)))))|^-) |o) 1^-) |^bit(^(^(^.)))^ 

+ 7l-/3Ph-(/3!i't(/3bit(/(r(j)))))b-> |1> |j> |/3'^"(/(r(j-)))) • (H) 

• Uncompute Q^'* by mapping 

|j) |6) \x) ^ \j) \b) \k) |(-x) mod 2™) 
and applying Qj'*. Now map 

|j) \b) \k) \x) ^ \j) \b) l(-fc) mod 2") la;) 

and again apply the copy operation. We can discard the additional quan- 
tum registers now, since they are in the state |0) |0) again. 
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A close look at the states in (|1U|) . and Hll|) reveals that wc have successfully 
simulated the mapping 

^phasc 

Since in general f3^^ o ^ id, our approximation introduced an error in the 
simulation. This error is bounded by 

sup|/5^"(/?"*(/(a;)))-/(x)| < sup |/?^'*(/3'^"(y)) - y| 

xeD yelOA] 

<sup{|x-y||/3^'*(a;)=/3^'*(y)}. (12) 

For the most obvious choice for 

: [0,l)->{0,...,2'"-l} 

a good /S'^'* is 

/3^'*:{0,...,2"-l}^[0,l] 

x^x- 2-" + 2-™-i 



(13) 



(14) 



for which the error from H12|l is less than 2 ^. 

As an example consider ^P^asc _ After some calculation we get 

Q phase _ ^phasG ii ^ n—m/2 

3 Phase queries cannot efficiently approximate 
bit queries 

Let us consider the (simple) evaluation problem 

with F = {f : {0} [0,1]}- It is obvious that this problem can be solved 
on a "classical" computer by just one evaluation of /, i.e. compj(S', i^) = 1. 
The precision e only depends on the accuracy of the query - the function call 
returning /(O). 

On a quantum computer that uses bit queries to evaluate / also 

comp^''^"(5,F) = 1, (16) 

at least, if we choose /S*"* as in (|13|l and if the number of qubits m is chosen so 
that 2™ > i. 

If we are only allowed to use phase queries however, we can prove a lower 
bound theorem. 

Theorem 3. Consider the evaluation problem For any quantum algorithm 

that is only allowed to access f by phase queries, 

comp«'f'"'"^(5,F) = r!(i), (17) 

if the encoding j^J'^"-'"^ is chosen as pp^'^^'^ = id. 
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To prove this theorem we will need lemma El and lemma |H1 thus we will 
postpone the proof. But theorem O allows us to easily prove the following 
corollary: 

Corollary 4. The approximation of an m-bit query Q^** to S < ^ by a phase 
query Q^'"^'"^ requires ri(2™) phase queries Q^j^"'^'^ when fjP^"-"^ — id. 

Proof. Suppose we could approximate an m-bit query by Uq = nq{m) phase 
queries and for all c > 0, mo € N there exists an m > mo with 

ng{m) < 02™. 

Fix such an m. Approximation means 



< s. 



We know that just one m-bit query is sufficient to solve the evaluation problem 
to accuracy e 2^'"+^ and certainty, so we could use the following algorithm: 

where 

UiQfUo 1^) 

is the algorithm that uses just one bit query. 

The probability that the output of this algorithm with an approximated 
query is incorrect would be less or equal than 2S, see (0). Thus we could solve 
the evaluation problem with 

comp^(S', F)<nq< c2™ = 2c^, 

which is a contradiction to theorem |3| □ 

With the result from section |21 that any phase query can be approximated 
by 2 bit queries we can conclude that bit queries are more powerful than phase 
queries. 

To prove theorem |3 let us introduce multivariate trigonometric polynomials 
(for an introduction see e.g. [S]). 

Definition 5. A mapping T : M" C is called an n-variate trigonometric 
polynomial, if it is of the form 



T{9o, . . . ,9n-i) = y^c,-e 



■i(n3.o9o + ---+"i,7i- 

■3^ 



with Cj G C, Uj^k G o,nd \ J\ < 00. We additionally define 

degr(6'o, . . . ,6'n-i) := max(|nj,o| + . . . + |nj,„_i|) 

as the degree ofT{6o, . . . ,6n-i)- 
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This definition allows us to define for every quantum algorithm, which de- 
pends on a phase query Q^^^'^, a corresponding multivariate polynomial. The 
idea that allows this is that every phase query as in Q can also be written as 

Of^^'=|j)|0) = cos0,|j)|O)+sin0,|j)|l) 
Qf^''\j)\l) = -sin^?,|j>|O>+cos0,|j>|l), (18) 
if 9j is chosen so that 



0j = arcsin ^^phasc(j(^Q-)))^ (ig) 

Now we can prove the following lemma, which closely follows the proof idea in 
Beals et al., 0. 

Lemma 6. For every quantum algorithm A making Uq queries Qj'"^*'^ there are 
trigonometric polynomials Tk{0o, ■ • ■ , 0n-i), 'iegTk{9a, . . . , 6'„„i) < Uq, so that 

A IV) = Un^Qf^^'Un^^i . . . U.Qf^^'Uo IV) 

= ^rfe(0o,...,0„-i)|fc) 

A; 

holds with 0j defined as in H9\) . 

Proof. The proof is by induction over the number of queries nq. 

Uq — : we have not made any queries therefore the state can be 

written as 

C/o|V) =^afe|fc) 

k 

with constants ak- 



Uq i-^ Uq + I : suppose 

C/n,gf . . . U^Qf^'^Uo IV) = ^ Tfc(0o, . . . , ft.-i) \k} 

k 

with degTfc(6'o, . . . , dn-i) < nq. We can split this state into 

^ rfc(0o, . . . , 0„-l) \k) = TrjM^ • ■ • . Sn-l) N) |j) |0) \l) 



-Y,T^M0o,■■■,e„.^)\^) \j)\l) 10 



(20) 



Let 9 = [da, ... , 0„_i). After applying Q?^'*^'^ to state (jSOJ, the result can be 
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written as 

J2 N> (cos^,T,,o/(e) \j) |0) + sin9,T,MQ) \j) |1» 10 

+ N> (-sin0,r,,i,(e) |j) |o) + cos0,r,,u(e) |j) |i)) \i) 

= (cos0,T,,o/(e) - s\ne,T,,um \i) \]) |0) |0 
+ {sme,T,MQ) + cos9,T,,nm) \t) \j) |1) \l) . 
The terms 

cos 9 jTijoi{9o, . . . , 6n-i) — sm6jTijii{6o, . . . , 6'„_i) 

and 

sin6jTijoi{9o, . . . , 6n~i) + cos6jTijii{6o, . . . , 6'„„i) 

again are multivariate trigonometric polynomials and of degree at most nq + I. 

Finally the mapping Un^+i is independent of /, thus not dependent on O, 
and just maps to a new linear combination of trigonometric polynomials, which 
again has degree at most Ug + 1. □ 

For (univariate) trigonometric polynomials we cite (see e.g. 0], 

Theorem 7 (Bernstein's Inequality). Let t{9) he a trigonometric polyno- 
mial. Then 

max \t'{9)\<degt{9) max \t{9)\. 

— n<9<7T —7T<9<7V 

The trigonometric polynomials Tk{9o, . . . , 9n~i) we found in lemma were 
related to / via equation 1)19(1 : 



9j = arcsin 



in^/3Ph-e(/(^(j))). 



What is the bound of Bernstein's inequality in dependence on arcsin ^? 

Lemma 8. // there exist an x ^ [0, 1], and a A with x + A (z [0, 1] so that for 
the trigonometric polynomial t{9) 



t(arcsinv^) > — and i ( arcsin VaT+A) < — 



and t{9) e [0, 1] for all 9 e R, then 



degt(^)>c(^^-+ V ^ j, (21) 

with m E {x, x + A} so that |m — i| is maximal and a constant c. 
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Proof. Applying the mean value theorem yields a ^ G [x, x + A] or ^ e [a; + A, a;] 
with 



|i (arcsin \/x) — t (arcsin V x + A)\ ~ arcsin yG: — arcsin \/ x + A 
and thus 

1 



— 2 arcsin \/x — arcsin y/x + A 

Now bound the degree oft{d) by Bernstein's inequality (lemma|7|) and use that 
me [0,1]: 



arcsin \^ x + A ~ arcsin ^/x 



<degt{0) max \t{e)\ < degt{e). (22) 

— 7Z<6<7V 



For the first bound we use that for all (fjtp E [0, 7r/2] 
-\(p- M < 1/21 sin^ If - sin^ VI, 
and therefore from (|22(l and 1)23(1 



(23) 



degt{e) > c. 



sm arcsm 



m Vx + A 



sm arcsm vx 




The second bound is proven by another application of the mean-value theo- 
rem. There exists a ^ e [x, x -I- A] or ^ e -I- A, x] with 

degi(6') > i (arcsin ^/iy{x + A-x) 

Note that for £_ £ [x , x + A] or S, e [x + A, x] 

(arcsin v^)' = 



is maximized by ^ = rn, with m € {x, x + A} so that |m — ^| is maximal. Thus 

, , \/m(l — m) 

degm > c^^^^ — '-. 

□ 

Note that the bound (|21|) closely resembles the bound of Nayak and Wu, [7] . 
With the help of lemma |H1 we are able to prove theorem O 

Proof of theorem\^ Let riq be the number of queries of an algorithm for the 
evaluation problem. By lemma trigonometric polynomials Tk{0) exist with 



(24) 



for e = arcsin A//3phase(/(r(j))), see and degTfc(6') < 
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The state (|24|) is measured and ip maps the resuh of the measurement to an 
element of our solution set [0, 1]. Let B be the set of all states so that 

B := {fc| |(p(fc) - 1/2| < e}. (25) 

The probability to measure a state \k) is given by \Tk{9)\'^ ■ Let 

I 

Then 

fees keB \ I J \ I ) 

which is a trigonometric polynomial again, degT(0) < 2nq, and T{d) G [0, 1]. 

Choose the input h.h^F with /i(0) = i, /2(0) = i - 2e, < e < i. If 
the algorithm is correct, 



T(arcsiny^/?phasc(j^(0))) > _ and r(arcsin y^/3phasc(j2(0))) < - 

(all query inputs are mapped to by r). 

Our lemma IHl tells us a bound for the degree of T(&): 

11 \/ mil — m) \ 
2n, > degT{0) > c U— + ^^-j^^ ^ j , 

where 

A = /?P^^^°(/i(0)) - /3P'^^^'=(/2(0)). (26) 

and m € {/3p^''"°(/i(0)), /3Ph^'''=(/2(0))} so that |m - \\ is maximal. 
We chose P^^"^" = id, so 



^(i_2e)(l-i + 2e)^ 



, ,1 v^2 - ^^n^ - 2 ^^^M / /I Vi""^^^ 

2n„ > c A / — 



2e 2e I I V 2e 2e 

There exists a constant c' so that for e small enough 

,1 

nq > c - 

holds, therefore 

comp^'P'^"''°(^,F) =17(i). 



□ 



4 Another proof for the hardness of evaluation 
with phase queries 

An anonymous referee suggested a shorter proof for theorem |21 The proof is 
based on the observation that for two inputs /i and /2 which are chosen as in 
the proof of theorem hJ2 & F = {f : {0} ^ [0, 1]}, 

/i(0) = i and /2(0) = i-2e, 
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the difference of Q^^^'^ and Q^^^'^ is bounded in the operator norm by 

|x),|||x)||=l 

= sup||(gP;^^^'=-QP^^^^^)(aoo|0) |0) +aoi |0) |1) || 

with ooo, ooi G C and |aooP + |aoiP = Ij since (5j|^^''° — Q^^^'' — for aU other 
inputs. 

Each of the terms can be bounded in the following way (here for aoo |0) |0)): 
11(0?^^'= - Q^'^aoo |0) |0) II 
= laool ||yT72|0> |0) + yT72|0> |1> - ^1/2 + 2e |0) |0) - v/l/2-2e|0) |1) 
I aoo I 



V2 
I aoo I 
a/2 



(1 - VT+ie) |0) + (1 - VT^e) |1) 



1 - 2\/l + 4e + 1 + 4e + 1 - 2^1 - 4e + 1 - 4e 



= |aoolV2 - a/TTI^- VT^e. 



If we develop + Ax into a Taylor series we get 

Vl + 4a; = 1 + 2a; - 2a;2 ± 0{\x\^) 



and thus 



||(gPf'^^^-QP;'^^'=)aoo|0>|0)|| 
< |aoo|v/2-2 + 4e2^0(|e|3) 



= laool ev/4TO(|e|) 
e |aoo|0(e). 

The bounds for the term aoi |0) |1) are analogous. 
Therefore we can conclude that 

||Qphase_Qphasc|| ^ 

Now we turn to our quantum algorithm, which is of the form described in 
equation |5J|: 

;7„gf^^'=[/„_i...c/iQf^^°c/o|v). 

We would like to bound the difference between the algorithm on input /i and 
on /2: 

f/„g^^^°t/„-i . . . U^Q}^^'Uo m - UnQt'^'Un-i . . . U^Qf^^'^'U^ II . (27) 
We will use the following inequality for unitary matrices A, B, C, and D: 

\\AB ~ CD\\ = \\AB -CB + CB- CD\\ 
< \\AB- CB\\ + lies - CDII < ||A - C|| ||B|1 + ||C|1 |1B - D|| (28) 
= ||A-C|| +||B-Z?|| 
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(see e.g. Nielsen and Chuang 8 ). Hence equation (|77|) is bounded by 



-jphasc ^phasc 



Recall the definition of B, the set of all states representing ^ from equation 
H25|l . Our algorithm for the evaluation problem must ensure that for /i the 
probability to measure a state from B is greater than | and for /2 less than j. 
Let p/j^ be the probability to measure a state from B on input /i and p/^ be 
the probability on input /2. We want to have 



\Pfi -P/2I > 2' 

We use the fact from equation Q that for A f\ : 
UnQ^f^^'^ ■ ■ - Uq and a measurement projection P 



(29) 



Un 



^phase jj^ 



< 



(^1 {4^P{Aj, - A^,) + A\PAf, + {A\.^ - 4jPAf, - A^^PAj 

{^\ 4p(% - IV) I + 1 (VI (4 - Ai)PAf, m 
< \\PAf, m \\{Af,-Af,)m + \\{Aj,~Af,)m \\PAf, m 

<2||%-%||, 

where we used the Cauchy-Schwarz inequality. This bound can be found, e.g. 
in Nielsen and Chuang -H?:- 

We choose P as a projection on the subspace spani? and get 



\pf, -pfJ<2\\Af,-Af, 
<2n gP)^""''- 



-jphase 

-/2 



(30) 



e 2nO{€) 

Combining equations (|29|) and Ij^UI) now yields 

n e f2 ^ 
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